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Abstract. Recently, an increasing interest in astrophysical as well as laboratory plasmas has been manifested in reference to 
the existence of relativistic flows, related in turn to the production of intense electric fields in magnetized systems [1]. Such 
phenomena require their description in the framework of a consistent relativistic kinetic theory, rather than on relativistic 
MHD equations, subject to specific closure conditions. The purpose of this work is to apply the relativistic single-particle 
guiding-center theory developed by Beklemishev and Tessarotto [2], including the nonlinear treatment of small-wavelength 
EM perturbations which may naturally arise in such systems [3]. As a result, a closed set of relativistic gyrokinetic equations, 
consisting of the coUisionless relativistic kinetic equation, expressed in hybrid gyrokinetic variables, and the averaged 
Maxwell's equations, is derived for an arbitrary four-dimensional coordinate system. 



A basic prerequisite for the formulation of a consistent relativistic kinetic for strongly magnetized plasmas in as- 
trophysical problems, is the formulation of single-particle dynamics in the context of a relativistic, fully covariant, 
formulation of gyrokinetic theory [2, 3, 4]. As is well known, this regards the so-called "gyrokinetic problem", i.e., the 
description of the dynamics of a charged particle in the presence of suitably "intense" electromagnetic (EM) fields re- 
alized by means of appropriate perturbative expansions for its equations of motion. The expansions, usually performed 
with respect to the ratio e = ri/L « 1, where L and ri are respectively a characteristic scale length of the EM fields 
and the velocity-dependent particle Larmor radius /"l = ' "^^^^ ~ Ik Larmor frequency and w the orthogonal 
component of a suitable particle velocity. The goal of gyrokinetic theory is to construct with prescribed accuracy in 
£ the so called "gyrokinetic" or "guiding center variables", by means of an appropriate "gyrokinetic" transformation, 
such that the equations of motion result independent of the gyrophase ^, ^ being the angle of fast gyration, which 
characterizes the motion of a charged particle subject to the presence of a strong magnetic field. In non-relativistic 
theory the gyrokinetic transformation can be constructed by means of a perturbative expansion of the form: 



INTRODUCTION 



z 



z = 



z'o + ez'i +eV2 + .. 
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which in terms of the Newtonian state x = (r, v) reads: 



r = r' + ep' + e^r2(z',f,e), 
V = M'b' + w' + V' + evi(z',f,e), 



(2) 
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where ep' is the Larmor radius. 
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while V is the electric drift velocity 



and the phyrophase gyrophase (j)' is defined: 
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In the past several methods have been devised to construct hybrid gyrokinetic variables. These include perturbative 
theories based, respectively, on non-canonical methods (see for example [5]), canonical perturbation theory [6, 7], 
canonical and non-canonical Lie-transform approaches [8, 9], as well as Lagrangian non-canonical formulations which 
make use of the so-called hybrid Hamilton variational principle [10, 4, 2]. 
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FIGURE 1. Guiding center and Larmor radius in non-relativistic theory. Here [e'j, 6^,63 = b'] denotes a right-handed basis of 
unit vectors. 



RELATIVISTIC GYROKINETIC THEORY: MOTIVATIONS 

For a certain class of problems in plasma-physics and astrophysics, existing limitations of the standard gyrokinetic 
theory [16, 17, 18, 19] make its use difficult or impossible. In particular, this involves the description of experiments 
in which the electric field may become comparable in strength to the magnetic field (so that the drift velocity becomes 
relativistic), and the study of relativistic plasma fiows in gravitational fields, which are observed or assumed to 
exist in accretion disks and related plasma jets around neutron stars, black holes, and active galactic nuclei. The 
finite Larmor radius effects and the influence of short wavelength electromagnetic perturbations are also expected 
to play a fundamental role in particle dynamics. In many respects, previous relativistic theory results inadeguate for 
such a task. In fact, some of mearher treatments consider the particle velocity as relativistic while its drift velocity 
is not [11, 12, 13, 14, 15]. This deficiency has been pointed out by Pozzo and Tessarotto [4], who developed a 
special-relativistic theory which includes the effect of relativistic drifts. However, the self-consistent inclusion of 
the gravitational fields, a prerequisite to make the theory suitable for astrophysical applications, as well the treatment 
of nonlinear EM perturbations of relativistic plasmas, requires a relativistic covariant formulation. This has been 
investigated by Beklemishev and Tessarotto [2, 3]. In this theory the appropriate relativistic gyrokinetic theory has 
been carried out through second order in the expansion parameter, including wave fields, based on a Lagrangian 
approach making use of the hybrid Hamilton variational principle. In such a case the variational functional for a 
charged point particle with the rest-mass ma and charge qa in prescribed fields can be written: 

S= J QtiAx^ = j {qA^,{x'')+u^,)Ax^', (7) 

where q = qa/maC^, and variations of occur on the seven-dimensional hypersurface UnU^ = 1, being the four- 
velocity = ^ and the relevant tensor notations are standard. Thus, g^v denotes the metric tensor components, 
characterizing the coordinate system (and the underlying space-time structure) which provides the cormection between 



the CO- and countervariant components of four-vectors (henceforth referred to as 4-vectors) — guvA^, while the 
invariant interval ds is defined as 

ds^=g^ydx^'dx\ (8) 

where the Greek indices are assumed to go through /i,v = 0,...,3. 



THE RELATIVISTIC GYROKINETIC TRANSFORMATION 

The object of gyrokinetics is to introduce a new set of phase-space variables (called the "gyrokinetic variables") 
such that the variable describing the rotation angle along the Larmor orbit (i.e. the gyrophase <^) becomes ignorable. 
This happens, by definition, when the Lagrangian (or, more generally, the functional) is independent of Once an 
ignorable variable is found, the number of corresponding Euler equations is reduced by one, and the new variables 
allow simplified numerical calculations, as the motion is effectively integrated over the fast Larmor rotation. The one- 
to-one transformation from the original set of phase-space variables {xf',u^) to the gyrokinetic variables is called the 
"gyrokinetic transformation". In what follows, we use the Lagrangian perturbative approach to find those variables, 
which is equivalent (in broad terms), to the Lie-transform method, though more direct. 

First, we assume that the curvature radius of the space-time and the gradient lengths of the background electromag- 
netic fields are much larger than the Larmor radius characterizing the particle path. However, we allow for existence of 
wave-fields with sharp gradients [kpi ~ 0(1), including k\\pL ~ 0(1),] and rapidly varying in time [(OPl/c ~ 0(1)], 
while such fields are assumed suitably smaller in strength than the background field. (We stress that, unlike in conven- 
tional formulations of gyrokinetic theory, this type of ordering is required in a covariant theory due to the reference- 
frame dependence of the ordering assumptions involving space and time scale lengths of the perturbations.) For this 
purpose we introduce the ordering scheme following the notation of [16]: 

Qndxf' = {u^+q{^A^+Xan)}6xf', (9) 

where e and A are formal small parameters (they should be set to 1 in the final results) allowing distinction between 
the large-scale background field A^, and the wave-fields given by afi. We search for the gyrokinetic transformation 
(y') = (x'",0,ju,M||) <-> (x",m'^) in the form of an expansion in powers of e: 

x^=x'^ + £ev;(y), (10) 

.s=l 

where (j) is the ignorable phase variable (gyrophase), and M|| represent two other independent characteristics of 
velocity (to be defined later), x'^ is the 4-vector "guiding center" position, are arbitrary 4-vector functions of the 
new variables {y') to be determined. We require that Ts are purely oscillatory in (p, i.e., the ^-averages of Ts are zero, 
as a part of the.x'^- definition. Note that the above descriptions of the new variables (x'",<^,/I,M||) will acquire precise 
mathematical meaning only as a result of the search for the gyrokinetic transformation. 

This search consists in applying the expansion (10) to the fundamental 1-form (9) and imposing the requirement 
that it is independent of in each order. A convenient framework is provided by projecting all 4-vectors and 4-tensors 
along the directions of a suitable fundamental tetrad (t, /,/',/")• i-^-' ™ orthogonal basis of unit 4-vectors so that the 
last three are space- like, and 

V^e,;,^,T^lh'^r = l, (11) 

where y^—gei^i^y is the purely antisymmetric tensor. As a consequence the four-velocity can be represented in the 
form: 

= w (l'^cos(j) + I'^sin^) + u^, (12) 

which can be also regarded as the definition for the gyrophase (j): it is defined as an angle in the velocity-subspace, 
where we introduce the cylindrical coordinate system. This definition is covariant with respect to transformations of the 
space-time coordinate system, which may change the vector components, but not the vectors themselves. Furthermore, 
we assume that w and are independent of (p . Validity of this assumption is justified by existence of the solution (at 
least for a non-degenerate Faraday tensor). 

The ^-independent part of the 4-velocity u is not completely arbitrary, but satisfies certain restrictions following 
from the requirement u^u^^ = I for all (j): 

Uh = u\\Ih + UoT:h, (13) 



ul = l+w^ + u\\'^. (14) 

Any two of three scalar functions w,Uo or M|| can be considered independent characteristics of velocity, while the 
third can be expressed via (14). It is now straightforward to eUminate from 5G' terms oscillating in (j) by properly 
defining displacements Tj. This task can, in principle, be carried out systematically at any order in the relevant 
expansion parameters (in particular in e). Thus, to leading order in e to eliminate the gyrophase-dependent terms 
in the fundamental differential 1-form one must impose constraint: 

+ Kqa^ — qr^F^v = 0, (15) 

where y denotes the oscillating part of y, namely y = y — y, y — {y)(i, is the gyrophase-averaged part of y and F^v is 
the EM field tensor. If the above requirement (15)is satisfied, the gyrophase is ignorable and the hybrid variational 
principle in our approximation can be expressed as dS" = 0. As a result, the (^-independent functional S" becomes 

S" = J {(^A'i, + ^qan + u\iln + UoT:n^dx'^' + 'fid(l)Y (16) 



where fi is the relativistic wave-field-modified magnetic moment, accurate to order e' and = y 1 + + M|p. 

The equations of motion, expressed in terms of the relationships between differentials tangent to the particle orbit, 
can be obtained as Euler equations of the transformed variational principle [2, 3]. Using the ^-independent functional 
(16) in the variational principle 5S" = defines the particle trajectory in terms of the new gyrokinetic variables 
(x'^,/i,M||,<^). This set is non-canonical, but further transformations of variables (not involving (j)) also lead to (j)- 
independent functionals and can be used for this purpose. 

THE RELATIVISTIC GYROKINETIC VLASOV KINETIC EQUATION 

The single-particle distribution function can be written in general relativity either in the eight-dimensional phase space 
^{x^,Uv), lJ.,v = 0, ...,3, or in the seven-dimensional phase space /(x'^,Mv), where only 3 components of the 4- 
velocity Uy are independent, so that 



^{x^',Uv)=f{x'',Uv)5{^Ju^u^-l)e{u'^). (17) 

The 5— function here reflects the fact that u^u^ = 1 is the first integral of motion in the case of the eight-dimensional 
representation. 

The kinetic equation in both cases retains the same form and yields the colUsionless Vlasov kinetic equation, namely 

dxt^ \ as J duv 

although in the 7-dimensional case v = 1 . 2, 3 only, while is the dependent variable. Here {duy/ds) is a function 
of independent variables {x^,Uv) found as the right-hand side of the single-particle dynamics equations. The kinetic 
equation can be multiplied by dsln this way it can also be represented in the parametrization-independent form as 
follows: 

df „ df 

■^dx^' + ^duy = 0, (19) 
dxf^ duy 

where the differentials are tangent to the particle orbit. 

Due to general properties of variable transformations it is obvious that any non-degenerate transformation of the 
phase-space variables {x^ , Mv) ^ (y') will lead to the same form of the kinetic equation 

|^dy = o, (20) 



where the differentials are tangent to the particle orbit. In particular, this property is useful for transformation to the 
gyrokinetic variables. 



Let (y) = Mil), then the kinetic equation becomes 



^fdx'^' + ^dm + ^dli + ^d(l) = 0. (21) 



By definition of the gyrokinetic variables the dynamic equations should be independent of (j) , i.e., expressions for 
{dx'^/d<j)) , (Jm||/J^) , {dp./d<j)) are independent of <j), while df/d<j) is periodic in <j). It follows that df/d<j) = 0, and, 
if ju is the integral of motion, d/i = 0, we get the kinetic equation expressed in the gyrokinetic variables as 

-dx'^" + -^dm = 0, (22) 



^x'^'■ du\\ 

which we shall call relativistic gyrokinetic Vlasov kinetic equation. Here the coefficients dx'^ and dM|| must be 
determined from the equations of motion in the gyrokinetic variables. 



THE MAXWELL'S EQUATIONS 

Finally we point out another important feature of the present formulation of the gyrokinetic theory. Namely, the 
Jacobian of the transformation is simple enough to allow exphcit integration in the gyrophase, needed for evaluation of 
the charge and current densities. The general form of the Maxwell's equations in presence of an arbitrary gravitational 
field is well known[20]. The first pair of equations can be written as 

ei^^^y^=0^ (23) 



while the second as 



where 



1 ^ (^^^=-1^/, (24) 



r-gdx"^ ^ ' c 

d^w 



/=c£^a/'«''/a(x,u)5(v4^-l)^ (25) 

is the current density, expressed via the distribution function of particle species a, and the signed particle charge qa- 
The 5-function under the integral allows to make partial integration, for example over mq, and arrive at a more widely 
used form 

o/ r^— ,\ .4 dMidM2dM3 

o V" "v-1 du^ . (26) 



However, in the gyrokinetic transformation the four-velocity is expressed via Eqs. (12)-(14) as 

u^i =w[l'^cos^ + l'^sm^) +u\\l^i+UoT^i, (27) 
so that d^ =wdwd^du\\duo [the sign is positive due to Eq.(ll)], while the partial integration over dug leads to 

5 (yu^'uv - l) ^A^^ ^^^y^W ^ (28) 



where Uo = y 1 + + My^. As a result, the expression for components of the current density can be rewritten as 

/=c£^a /(w(/'''cos0+/"^sin(|))+M||/''+M<,T'')/a(x,u)^:^^^^ff^. (29) 

Further, the distribution function fa is expressed as the function of the gyrokinetic variables 

fa=fa{x'^,%u{), (30) 



and it is necessary to transform it back to particle coordinates before integrating, as in Eq.(29) the particle position 
X, rather than its gyrocenter position x', is kept constant while integrating over the particle velocity. This makes it 
convenient to rewrite Equation (24) as 

±- {^gF^'^) = -^/^= eM(x), (31) 
where the right-hand side is also evaluated at x. Then 

Q^'{^^) = -AnY,qa I [w{l'^'cos^ + l"^'sm^)+u\f+UoT^]fa (x-£e'r; ) (32) 
a J V ,=1 / 



CONCLUSION 

A closed set of relativistic gyrokinetic equations, consisting of the collisionless gyrokinetic equation and the averaged 
Maxwell's equations, is derived for an arbitrary four-dimensional coordinate system. 

In several respects the theory here developed represents a significant improvement with respect to kinetic equations 

derived by other authors. The present covariant kinetic theory adopts a set of hybrid gyrokinetic variables, two of which 
include the Lorentz-invariant magnetic moment and gyrophase angle. The theory, allows E /B ^ 0(1) and therefore 
permits relativistic drifts (Vd ^ c) an moreover takes into account nonlinear effects of the EM wave-fields. Moreover, 
since the gyrokinetic transformation is obtained to the second order in terms of the ratio of the Larmor radius to 
the inhomogeneity scale, the theory can be applied also to the investigation of finite-Larmor radius effects. Another 
interesting aspect is that in the present theory the wave field is no longer limited in frequency and the wavelength, i.e., 
(0/Q.c ~ 0(1), A:||Pl ~ 0(1), so that the class of admissible waves is broader than the usual "drift- Alfven perturbations" 
and can include the magneto-sonic waves, for example. 
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